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We apply the functional renormalization group (FRG) method to calculate the conductance of a quantum dot 
in the Kondo regime. Starting from the exact FRG equations in Keldysh formulation for the Kondo exchange 
Hamiltonian in pseudo-fermion (pf) representation, we solve the coupled equations for the pf self energy 
and the coupling function, neglecting three-particle and higher correlation functions. The conductance G as 
a function of temperature T and bias voltage V is calculated using a renormalized Golden Rule expression. 
The limiting behavior at T and/or V 3> Tk {Tk '■ Kondo temperature) agrees with known results. The 
difficulties when approaching strong coupling are analyzed and improvements are suggested. 



1 Introduction 

The Kondo effect is a ubiquitous phenomenon in the area of electron transport in nanostructures 0] 12 El 
|4][5]|6). F° r i ts existence a small system with internal degrees of freedom (e.g. a quantum dot with spin), 
coupled to a sufficiently large conduction electron system (the leads) is required. If certain symmetries are 
satisfied, the localized spin forms a resonance state with the surrounding conduction electron spins at the 
Fermi level, of energy width Tk , the Kondo temperature (here and in the following we use units where 
Boltzmann's constant fee = 1 and Planck's constant h = 1 ) . A finite bias voltage V applied across the 
quantum dot is known to destroy the Kondo resonance, provided eV 3> Tk- It has been argued that the 
limitation in the coherence time brought about by inelastic processes in the leads causes the suppression 
of the Kondo effect at finite V., since it cuts off the renormalization group flow of the coupling constant 
[ 8 , 7 1 . At sufficiently large V this cutoff is still within the perturbative regime. In that case the perturbative 
RG treatment ("poor man's scaling", J9) ) generalized to nonequilibrium, proposed in ifTol . provides a 
controlled and complete description of all physical quantities. Indeed the latter theory has been used to 
decribe experimental data for various situations, each time with excellent success ifTTl [121 [131 [141 . While 
the effect of decoherence on the RG flow of the coupling constant has been included in iflOl QT) in a 
phenomenological way, it is of interest to derive the effect of the decoherence rate in a systematic fashion. 
The way in which the decoherence rate cuts off the logarithmic singularities in perturbation theory has been 
analyzed in 11301 . where it has been found that the decoherence appears in the form of physically accessible 
rates, namely the longitudinal and transverse spin relaxation rates. In a completely new way the effect of 
decoherence on the Kondo effect has been analyzed within the flow equation method 11151 . with results in 
agreement with and thus confirming the work by Paaske et al. lfT0ll30l . The latter method is, however, not 
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easily applicable to the calculation of observable quantities like the conductance. It is therefore worthwhile 
to work out a systematic and consistent derivation of the coupled RG-equations for the coupling function 
and the decoherence rate, even in the weak coupling regime. 

On the other hand, it has not been possible , so far, to calculate the influence of a finite bias voltage in the 
strong coupling regime. There are several methods available describing properties of the Kondo system at 
strong coupling in equilibrium (see [ 16 1), but a general and practically applicable extension of these meth- 
ods to nonequilibrium has not been found yet. A recently proposed "Conserving T-matrix approximation 
", based on self-consistent renormalized perturbation theory including all two-particle channels, which has 
been shown to recover correctly the Fermi liquid behavior at low energies and the crossover to the local 
moment regime 1171 [T8l . may in principle be generalized to nonequilibrium, but is at present too costly in 
numerical evaluation effort. 

In this paper we present results for the weak-coupling regime of the Kondo model in and out of equi- 
librium obtained by the functional renormalization group (FRG) method [19 20] [21] [22]. The FRG is a 
powerful tool when dealing with the diversity of energy scales and the interplay of correlations in an inter- 
acting fermionic systems such as the Kondo model where the microscopic energy scale band width D and 
the dynamically generated energy scale Kondo temperature Tk differ by orders of magnitude. The FRG 
uses an exact mapping of the initial model to models with reduced band width. The mapping generates 
an infinite hierarchy of flow equations, describing the gradual evolution from a microscopic model Hamil- 
tonian to the effective action as a function of a continuously decreasing energy cutoff. In practice, one 
modifies the bare propagator by cutting off frequencies smaller then some infrared cutoff A. As this cutoff 
is progressively lowered more and more low-energy degrees of freedom are included, until the original 
model is recovered for A — * 0. The conventionally used approximation consists of a suitable truncation of 
the hierarchy of equations. We find that the truncation neglecting three particle and higher order correlation 
functions, and keeping only the self energy together with the four-point vertex function captures the weak 
coupling regime (T, V 3> Tk) completely. In this paper we show in detail how the transport properties of 
a Kondo quantum dot in the weak-coupling regime may be calculated in a systematic and complete way. 

Compared to previous applications of the FRG method to Kondo physics ( 11231 and ref. therein) our 
approach differs in at least two ways: (1) the present work addresses the Kondo model, rather than the 
Anderson model (compare, however, the recent work 11241 in which the dynamical spin response of the 
Kondo model has been considered, using a real-time RG method [33]); (2) in contrast to previous FRG 
works we include the effect of the decoherence rate. 

2 Model and FRG formalism 

We consider the model of a quantum dot described by a spin | operator S coupled by exchange interaction 
and tunneling to the conduction electron spins as described by the Hamiltonian 

H = Y^^ ~ ^)c\ aiJ c kaa + J S 'W 

kaa a.a 1 — 1,2 

Here a, a' = 1, 2 labels the leads assumed to be in equilibrium at chemical potentials = ±eV/2 and 
a =t, I are spin indices. The conduction electron spin density operators at the quantum dot are defined as 
s QQ ' = \ c \aa Ta , c " c k'aia' > w ^ c feacr an electron creation operator in momentum state k, lead a 

k,k' ,cr,a"' 

and spin state a and r the vector of Pauli matrices. The s QQ/ are tunneling operators for a ^ a'. We use 

a pseudo-fermion representation for the local spin S =| l\ T ~a^~ii > with projection onto the physical 

7,7/ 

Hilbert space by adding a term \[%2 f- f f- f ~ 1] t0 me Hamiltonian and taking the limit A — > 00 (cf. 112511261 ). 

a 

The FRG equations for the self energy and the one-particle irreducible Green's functions are derived by 
mapping the generating functional for the connected Green's functions in Keldysh space (see ll20l |2"T1 l22l 
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for an introduction into the FRG method) 

W[r], x) = log / D[xj), </>] exp[iS] exp[(-0, 77) + X ) + h.c_ 



with the action 

S = (rjj, (G )- V) + (6 *u~V) ~ 0] 

where Go and Fo are the free propagator Keldysh matrices for the lead electrons and pseudo fermions, 
respectively, and V[^>, cf)] denotes the exchange interaction, onto a generating functional in which certain 
high energy degrees of freedom have been removed at the expense of changing the interaction part accord- 
ingly. We use standard notation for the functional integrals (see e.g. IF271D . In particular all time integrals 
are extended over the Schwinger double time contour. 

Since the interaction is local in space we may integrate out spatial degrees of freedom. Thus, for station- 
ary situations, the propagators depend only on a single frequency. We choose to remove the high energy 
components of the lead electrons, by introducing a sharp and symmetric — with respect to fi a — cutoff A , 
expressed by the step function 0(|w — fi a \ — A), according to 

G° Q A M = G»e(|c- Ma |-A) 

and resulting in a A-dependent generating functional Z A . It is convenient to perform a Legendre transfor- 
mation to the generating functional W A for the one-particle irreducible (1PI) n-particle vertex functions 
I A . An expansion of the functional T A in the fields I A leads to the flow equations for the I A , and in 
particular for the self energy £ (we use £ instead of 



and for the two-particle vertex function 1 = 1 



2 ■ 

o raac.a'a'd/ i i\ 1 \ ^ T aac,a<7c' / ^0,d'c / / \ -j-ib' a' { , \ -rCicrd'oL <j' d 

d^ya^'b (We.W/iWe.Wj) = -5- ^ G n (e)F (u e + UJf - eJJ^ Jj, 

. T aad' ,a' cr' d,-~,0,d' c' r \ jpb' a' / , , , ,1 , \ T acrc,aac' 
+V,7a' G S \ e > h \ U S ~ UJ e+ e ) 1 jb',-y'b ( 2 > 

For the sake of simplicity we occasionally omit the frequency arguments of the vertex function. Keldysh 
indices are denoted by the latin indices a,b,c, d. At this stage we have already neglected the term in equ. 
(0 which involves the three particle vertex function ^3 . This term would obey the next RG equation in the 
infinite hierarchy of equations for the /„ . The principal approximation of our approach is the neglect of 
^3 and all higher vertex functions. In Fig. [T]the flow equations (Q~|) and (ffl are depicted graphically. The 
boundary conditions at the initial value A = Ao = Dq are given by 



= (3) 

T OLOC.OLo'd( f f\ -Jcxa.' -> -> . c d / a\ 

4a,a'b K,W/;W e ,W/) = I— — WT 77 'A ab , (4) 



independent of frequency. Here, A^ denotes the bare Keldysh vertex structure, which is given by 

Kb = (-l) a+1 5abS cd or k a d b = \(5 ab Tl d + Tl b 5 cd ) (5) 

in the contour ordered or the Keldysh rotated representation, respectively. Furthermore we use either a 
representation of the Keldysh matrices in the usual Schwinger contour ordered representation (T, T, <, >) 
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n 



Fig. 1 Graphical representation of the fRG equation for the self energy E and the two-particle vertex function 7 as 
given by equ. (TJ and (2). The dot on the l.h.s. represent the derivate w.r.t. A. Solid lines with a stroke display the lead 
electrons propagator single scale propagator S A = G°' A ■£^(G°' A )~ 1 G°' A , whereas the double dashed lines represent 
the dressed pseudo-fermion progagator. 



or in terms of retarded, advanced and Keldy sh components (R, A, K), respectively. In both representations 
the matrix Green's functions are given by ll28ll 

G T G<\ . _ (G R G K 



G= {g> G f ) ° r G= {0 G A ) (6) 

Corresponding representations are used for the pseudo-fermion self energies. The bare conduction electron 
local Green's functions at the quantum dot in both leads are given by 



M = -i*N e(D - M) = [G"- (w)]* 

G° a K (u) = -2«tanh {^^j N @(D - \w\) 

G° a <(uj) = i2irf a (w)N Q(D - \u\) 
G%>(u>) = -*27r(l - f a (u))N Q(D - |w|) 

where iVo = 1/2Dq is the density of states, f a (uj) = l/(exp((w — /j, a )/T) + 1) is the Fermi function. 
Since we do not include a Zeeman term for the conduction electrons, the lead electron Green functions do 
not depend on spin and hence the spin indices have been suppressed. The dressed pseudo-fermion Green's 
functions read 



F« x {u) = -iA 7 (w)(l - 2n 7)A ) 

F> x {u) = -i{l-n ltX )A^u) 
with the spectral function 



[u + a§-R y (u>)] 2 + [r a (u>)/2}2 
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Here the imaginary part and the real part of the pseudo-fermion self energy are given by Y = i(E R — 
T, A ) or R = Rc(S fl ) , respectively. The pseudo-fermion occupation numbers are denoted by n lt \ . 
Note that n 1: \ — > in the limit A — > oo (compare the calculation of Y in ll26l ). In the presence of 
a magnetic field, however, it is necessary to keep the n^ t \ , as they describe the local spin polarization 

m = lim A ^oo[e A/T (n T! A - n^x)]. 

3 FRG method applied to the Kondo model 

3 . 1 Vertex function structure 

Within the perturbative RG method as formulated in iflOl [TT1 only a single component of / in Keldysh 
space has been kept. We now consider the structure of the vertex function / in Keldysh space in more 
detail. Starting point for the following considerations is equ. d2J in the Keldysh rotated form. The leading 
logarithmic terms in the FRG equation for the vertex function originate from the product of the Keldysh 
component of the lead electron Green function G K with the real part of the retarded or advanced pseudo 
fermion propagator F R ' A . Keeping only these terms in the sum over Keldysh indices on the right hand side 
of equ. (O and neglecting the A dependence in the energy arguments of the vertex function this equation 
may be written in a short hand notation (keeping only Keldysh indices) as 



rid T c2 



/ j '^R,A I aa/ I a. 



(8) 



Here n and C denote the kernels in the Peierls and the Cooper channel, respectively, and the leading 
components derive from its imaginary parts. For completeness, however, we state both the imaginary as 
well as the real parts. To this end we decompose the individual channels as 



with 



and 



Cr a — iC r + iCi 



J R,A 

n% A = ±n r 



lYi,, 



C, 



tanh(^)(w/ +w e -e-/j a + 7§ ) 



2irN 

IT 

27T7Vn 



E 

e=±A {bJf + U e — € — (la 



7f) 2 



2 



/ 



e + /i a 



iff 



tanh(2^) Y 



c=±A (uif + uj e — e - Ma + 7 -j) 2 + [ 
y-v tanh(^)r 7 

e=±A (w'f - I 



/'-o 



7#) 2 + ^ 



(9) 
(10) 



(ID 



(12) 



(13) 



(14) 



To simplify matters we have already assumed an energy independent imaginary part of the pseudo fermion 
self energy T 7 and furthermore dropped the real part of the self energy. Again we emphasize that C, , oc 
j- and hence produce the leading logarithmic components, whereas C r , H r oc give rise to sub-leading 
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corrections. In the following we neglect contributions from C r and II r and focus on the leading compo- 
nents only. Hence, we are led to the equation 



a/=l,2 a/=l,2 



Id 

aal ± alb 



(15) 



To examine the behaviour under the FRG we perform a Picard-Lindelof like iteration, i.e. we start with 
the bare vertex structure I£g — IA.°f = ^I(5 a bT^ d + T^ b 5 c d) on the right hand side of equ. (1151 1. By 
performing the sum over Keldysh indices we see that in both the Cooper and Peierls channels this structure 
is reproduced. Hence, in leading order there is only one invariant amplitude, of the structure of the bare 
vertex 

Ia d b~IK d b- (16) 

confirming the treatment in. iflOl [TP . To recover the perturbative RG scheme of iflOl ITT1 from the FRG 
equation (O the following additional steps have to be performed: As a first step the imaginary part of the 
pseudo fermion self energy is neglected in the denominator of equ. (TTTT) and (fT2l) . Moreover the energies 
of the pseudo fermions are put on-shell, i.e. Uf — —jB/2. Then the coupling function depends only on 
a single frequency, the frequency of the incoming electron.. Finally the cutoff dependences in the Cooper 
and Peierls channels (equ. (fTTl i and ( fT2b ) are approximated by a window function 



\ - sign(e) _ A_ 



^e-Auj A 2 - (Aw) 2 

e=±A v ' 



26(A-|Aw|) (17) 



as described in more detail in IflOl ITTI . Here Aw depends on V, B and the incoming and outgoing fre- 
quencies. Finally, in order to incorporate decoherence effects, which serve as a cutoff for the RG flow, a 
relaxation rate Y , which has contributions from both, self energy and vertex corrections, is incorporated 
in equ. ( fTTT i by replacing 0(A — |Aw|) by 9 (A — y/ (Aw) 2 + T 2 ), where Y in turn is calculated by a 
second order golden rule expression with renormalized vertices. These two equations (i.e. for the coupling 
function and the decoherence rate Y) are then solved iteratively, until convergence is reached. 

In the following we will show how the decoherence rate (to be more precise the imaginary part of the 
pseudo fermion self energy, which is part of the decoherence rate) arises in a natural way. To this end we 
first simplify the coupled equations (TJ and (fJI in an appropriate way and then solve these equations for 
different parameter regimes. 



3.2 Approximated flow equations 

In order to simplify the coupled differential equations for the imaginary part of the pseudo fermion self 
energy Y and the effective interaction I we first neglect the frequency dependence of the leading component 
of / , 



o t ol(T\ol' a' 1 \ ^ I" T ct(T,ctcr rOifrio. ' <j' s~ia. i a \ i 7-acr, aa T-acriacr-i-rQ /* \ 

d A / 7;7 , =-^2^[ I rn 7 7,7' C 7 (A)+/ 7; V ^ n?(A) 

CK7<T 



(18) 



where 



C|(A) = -2nN £ A* +7f) 



£ =±A (_ e - Mdi +7f)2 + (^y 
nf (A) = -2nN ^ tMlh(#)(£ + Ma+ ^ f) 

e=±A ( £ + Ma+ ^|)2 + (h. 



(19) 



(20) 
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The RG equation for the imaginary part of the pseudo fermion self energy T = i(Y> R — J^ A ) requires 
knowledge of the vertex function I"^^" 7 (e, 0; e, 0) at e = ±A. By simply replacing this vertex function 
by an energy independent one of the form ( fT~6b the imaginary part of the pseudo fermion self energy is seen 
to vanish identically. Here it is necessary to take into account the subleading Keldysh components of / . 
To this end we use the leading components on the right hand side of the equation for I(e, 0; e, 0) , equ. (0, 
to generate the subleading Keldysh components, which we denote by 1(e) : 



■ ;i jaac-.aad / \ ^ 



acr.acr j-Q(7:acr 



TLXU 

7;7 



T 

717 



a cc' \d dfic d ot f k _\ I jolo .ol<7 jolg\ol<7 a d d \ cc jtc' d' a / a \ 
A aa' y V6°a'6'7 \ l \ e )^~ 1 ra J 7 ; 7 A aa' A 6'b U a'6'7 l A > 6 ) 



(21) 



where the Einstein summation convention is used. The kernels in the Cooper and Peierls channel read 

Cii* (A, e) = Gi c ' (e)if °' (e - e) (22) 

£ = ±A 

IT^4(A, e) = J2 Gt c ' (e)<°' (e + i) . 



(23) 



e=±A 



With the aid of (f2Tb the equation for the imaginary part of the pseudo fermion self energy may be cast in 
the form 



(24) 



e=±A 



Hence equ. ([T8T l. (ETb and (124-b provide the desired system of coupled differential equations for the vertex 
function and the imaginary part of the pseudo fermion self energy. For the subsequent analysis we are 
going to use dimensionless couplings g = NqI. Additionally, the spin dependence of the vertex function 
is parameterized as g. 
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Ty~/i and, for the sake of simplicity we restrict ourselves to the 
case of symmetric couplings g aa i = g. Furthermore, in order to establish a more direct connection to 
Poor man's RG schemes lfT0lfTTl l9ll and to prevent a factor of i from pervading the formulas, we employ a 
redefinition of g, i.e. g — > i8g, which results in the new boundary condition g = go = No J at Ao = Do- 

3.3 Weak coupling and linear response at finite temperature 

To begin with a quantitative analysis we focus on the weak coupling regime B, V — > and T/Tk ^> 1, 
where the Kondo temperature is defined by Tjc = Do exp(— l/(2go))- In this case the Cooper and Peierls 
kernels $1% and d20l i simplify to 



C(A) = 4tt7Vo 



Atanh(^) 



-n(A) 



A 2 + (r 7 /2r 

Performing the internal spin sums in equ. ( fTST l , one finds 



d\g = -2.g 2 tanh ( 



A 



2T J A 2 + T 2 ' 

Similarly, in the absence of a magnetic field the RG equations for V and g take the simpler form 
5A5(e) = 3. 9 2 £ /(-?) 



€ = ±A 



(r/2)' 



d A T = £ f(e)g(e) 



(25) 

(26) 

(27) 
(28) 



e=±A 
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2.5 




1 10 100 



T/T K 

Fig. 2 Leading component of the coupling func- 
tion g(0) as a function of T/Tk for values down 
to T/Tk » 1. Dots: results from FRG; solid line: 

9 = 21og(T/T K )- 




1 10 100 



T/T K 

Fig. 3 Imaginary part of the self energy F(0) 
versus T/Tk for values down to T/Tk ~ 1. 
Dots: numerical results from FRG; solid line: F = 

37rT (21og(T/Tjf)) ;i ' 



In this case the corresponding initial conditions read 

g Ao=D (e) = and T Ao=D = . (29) 

For the numerical treatment of these equations, however, we have to take a small but finite value for r A ° . 
It turns out that the solution of the FRG equations does not depend much on the choice of the initial value 
of r A ° as long as it is taken to be sufficiently small. 

We now solve the FRG equations for values of T/Tk 3> 1 and obtain results shown in Fig. [2] and Fig. [3] 
In Fig. |2]the leading Keldysh component of the coupling at the Fermi energy, g(u> — 0), is shown at zero 
bias as a function of T/Tk- At temperatures T ^> Tk the numerical result (dots) agrees very well with the 
result of renormalized perturbation theory, g — 1/2\tl(T/Tk) ■ It is seen that the effect of spin relaxation 
leads to a reduction of the growth of g(u> = 0) as compared to the weak coupling result, as expected, 
but does not cause g(ui — 0) to take a temperature independent value of order unity at T <C Tk ■ The 
relaxation rate V , shown in Fig. [3] also agrees well with the renormalized perturbation theory expression, 
T = 37rT / '[21n(T '/T K )] 2 at T > T K . At lower temperatures T is seen to fall well below the diverging 
perturbative result, as expected, but does not appear to approach a limiting constant value of order Tk at 
T <C Tk ■ 

The conductance G (in units of e 2 /irh ), in the absence of a magnetic field, may be approximately 
obtained from the golden rule expression 

G(T) = -g 2 {u = 0) (30) 

Fig. |4] shows a comparison of the linear conductance as obtained from equ. ( f30b (dots) and the weak 
coupling result (solid line) G = (37r/16)/[ln(T/TV)] 2 • Even though the rapid rise of G as T approaches 
Tk is somewhat slowed down by the effect of the relaxation rate T , the slowdown is too weak or happens 
too late so that G rises beyond the unitarity limit. 

It is somewhat surprising that the initial strong growth of g is not cut off more effectively by T . After all, 
at T > Tk the relaxation rate reaches values of T ~ 2>Tk so that it might be expected to be more important 
than the temperature in cutting off the RG flow. The reason for this different behavior is apparent from 
Fig. [5] which shows the behavior of V as a function of the RG cutoff scale A : A sizeable value of T is 
only generated at the scale A < T and not at Tk as one might have expected. Therefore, at T <§; Tk and 
A ~ Tk the relaxation width V is still exponentially small and does not serve to cut off the flow of g to 
large values, which explains why g, T,and G (at A = 0) appear to grow beyond the unitarity limit. We 
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emphasize again that the above results depend sensitively on the approximation of neglecting the energy 
dependence, both, of the coupling function g(u>) and the relaxation rate T(u>) (see below). 



3.4 Weak coupling at finite bias voltage 

As it is known iflOl [TT1 that in case of a finite bias voltage V the frequency dependence of the coupling 
function is crucial, in the following analysis we follow IfTOl [TP and maintain the frequency of the lead 
electron. We approximate the kernels in the Peierls and Cooper channels by using 

E tanh ^) ( £ - A^2 +V/ 2)2 -2tanh(A)e (A -V(A-) 2 + (r/2) 2 ) (31) 

where Ao; = u> — fi a . Hence, the equations for the leading component of the vertex function, g(u>) , the 
subleading vertex component g , and the imaginary part of the pseudo fermion self energy T take the form 

d A g{u) = -g 2 (u;)j ]T 8(A - y/(w - [i«] 2 + (W) (32) 

a 

W-^EftW p^, (33) 



2 

e=±A a 



where 



fa (e) = 77^Y-7 . Ha = ±eV/2. 

exp (— f 3 -) + 1 



Results of a numerical solution of these RG equations are shown in Fig. |6]and Fig. [7] In Fig. |6]we display 

g(uS) for various values of V/Tk- Fig. |7] shows the imaginary part of the pseudo fermion self energy V 

l 

\o S '(V/2T K ) 



(dots) with the approximate expression T = ^j-V-, — 'iiy/ot i ( s °lid li ne )> which is seen to agree very well 
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• V/Tk = 251.7 

• V/Tk = 111.8 

• V/Tk = 49.7 

• V/Tk = 22.1 




-200 -150 -100 -50 50 100 150 200 

uj/T k 

Fig. 6 Leading component of the coupling func- 
tion g (oj ) versus to / Tk for various values of V/Tk ■ 




Fig. 7 Imaginary part of the self energy V versus 
V/Tk- Dots: numerical results from FRG; solid 



line: Y = ^-V 



4 V (21og(V/2T K )p- 



in the regime V 3> Tk ■ 

At finite voltage the charge current I c is approximately given by 



3ir 2 e f du> f doji f de 



2H / 2tt / 2tt / 2tt 



■\ g (u J ,e;u;f,e + u~u;/)\ 2 [f L (Lu)-f R (u J )}A(e)A(e + ^^^) .. (35) 



In the limit of vanishing relaxation rate T , the above expression simplifies to 

Sire 



/,.= 



An 



dw\g^A^M 2 [fL{^)~ f R {w)]. 



(36) 



In Fig. [8] we show the conductance G = I c /V in units of Go as obtained from equ. d36*i l versus V/Tk 
at T = 0. Again, dots represent the numerical evaluation of the FRG equations, whereas the solid line 
represents an approximate analytical form G/Gq = ^rg~ \^(v/2T K ) • ^ ne obtains this analytical result by 
replacing g(u>) by g(oj = 0) in equ. d32b . However, because of the distinctive form of g(u>), which exhibits 
pronounced peaks at the resonance frequencies lo = ±V/2 the resulting G(V)/Gq is slightly enhanced 
compared to the analytical form. In Fig. |9]results on the conductance G/Gq versus V/Tk are presented 
for temperatures T/Tk = 50, 10, 5 , such that one remains in the perturbative regime even in the limit 
V — » 0. Dots are the FRG results, whereas the solid line represents G/Go = ^f-<? 2 with the approximate 
vertex function 

g(u = 0) w 1 . (37) 

i K 

The frequency dependence of g(u>) together with the integration over the difference of the Fermi functions 
in the left and right lead in equ. (f36t lead to a distinct discrepancy between the FRG result and the analytical 
form for V < T. In Fig. [I0]results on V versus V/Tk in the same temperature regime as in Fig. [9] are 
shown. As can be seen, in the temperature range 3K < T < 10K T saturates to a value of the order of Tk 
when the voltage falls below the corresponding temperature. 
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Fig. 8 Conductance G/Go 
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Fig. 9 Conductance in units of Go = 2e 2 jh as 
a function of V/Tk- Dots denote the result ob- 
tained by solving the fRG equations numerically. 
The solid line corresponds to the weak coupling 

/ V(V/2)2 +T2 



form G/Go = 3tt 2 1/16 log 2 




Fig. 10 Imaginary part of the self energy T as a 
function of V/Tk- Dots represent again the results 
from the fRG. 



3.5 Towards the strong coupling regime 

From the previous considerations, i.e. the analysis and the solution of the approximated set of FRG equa- 
tions [26] [27] and [28] it can be inferred that this set of equations in its present form is not appropriate to 
capture the physics in the strong coupling regime, i.e. in the regime T, V <C Tk- There are two major 
deficiencies in the approximate RG equations: (1) the unitarity limit on the conductance, G ^ Go is not 
respected; (2) the spin relaxation rate T is found to rise with decreasing temperature/voltage in the strong 
coupling regime, whereas Y is expected to saturate to a value of order Tk ■ There might be several reasons 
responsible for this failure. 

First, the effect of higher order correlation functions is neglected in the approximation scheme discussed 
so far. The additional terms in the RG equation for the coupling function g(u> = 0) generated thereby are 
similar, but not identical, to the higher order loop corrections to the poor man's scaling approach, and may 
be expected to slow down the growth of the coupling near the energy scale Tk ■ A quantitative control 
of these terms within FRG appears to be out of reach. Nonetheless, in the following we shall explore the 
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consequences of a model assumption including higher order terms. To this end we add terms of third and 
fourth order to the r.h.s. of equ. d26l i 



d A g = -2 5 2 (1 -ag- (3g 2 ) tanh ( ^\ , (38) 

with coefficients a, (3 to be determined by requiring that the unitarity limit be satisfied and that the resulting 
conductance fits the corresponding NRG result |29| best. We find that taking a = and putting j3 = g 2 = 
4/ (37r) , in equ. (|38l leads to the best agreement with the NRG data on the linear response conductance 
11291 . Note that the additional terms are not interpreted as two-loop and three-loop contributions to the 
(3— function. Rather, they are thought to represent the zero of the (3— function at the strong coupling fixed 
point. A complete discussion of the (3— function at equilibrium for the anisotropic Kondo model in the 
strong coupling regime has been given in DTI , on the basis of the flow equation method. The simple 
ansatz used here cannot be compared with the entirely different approach of Oil . The point we want 
to make here is that a simple adjustment of the f3— function for the coupling constant is not sufficient to 
capture the physics of the relaxation rate (see below). 

In the following we keep the frequency dependence of T(uj) . Within the pseudo fermion projection 
scheme using the limit A — > oo the self energy acquires a highly asymmetric form, which causes problems 
in the numerical solution. We therefore choose an alternative projection with A = 0, which requires to 
introduce normalization factors accounting for the difference in the impurity partition function with and 
without projection. In the calculation of r(w) in the absence of a magnetic field these factors are not 
necessary. This choice has the advantage of particle-hole symmetry, i.e. T(uj) = r(— u>). 

The FRG equations (cf. H~8| [2T1 and[24b then take the following form 

»A* = -9 2 [1 - {9/9uf] A2 + A p2(A) {^r) (35) 

d A I» = ^^[/ a ( f )j>( e ,,) + (l-/ a ( f ))s<( e ,,)] (40) 



2 

e=±A a 
e=~±A r - 



e =±A a 

Fig. [TT| shows a comparison of the linear conductance from FRG (dots) with the NRG data l29l . 
For the purpose of this comparison we adopted the definition of a Kondo temperature T£ used in [29 1 
: G(TJ-) = Go/2. This differs from our previous definition of Tk = Dq exp(— 1/(2<7o)) by a factor 
Tk/Tk ~ 2.76. We find excellent agreement between the FRG and NRG results. In Fig. Q~2] we show the 
imaginary part of the pseudo fermion self energy T(u> = 0) as a function of T/T^. The inset shows the 
region T <C T^, where T again decreases almost linearly with decreasing temperature, down to a similar 
residual value of T/T^ « 0.002 as found from the weak coupling FRG equation above. 

We now turn to the nonequilibrium situation: In Fig. [I3]we show the conductance as a function of the 
voltage. The lines indicate the half width at half maximum which amounts to AV — 1.93T^-. Fig. [14] 
shows r(o; = 0) as a function of V/Tg. In the inset T is seen to drop linearly with decreasing voltage at 
V < T£, down to a residual value of Y/T^ w 0.002. 

To summarize, adding a fourth order term with adjusted coefficient to the /3-function of the RG equation 
for the coupling g(u> = 0) leads to good agreement with the exactly known conductance in the limit V — > 0. 
The relaxation rate T(u> — 0), however, comes out three orders of magnitude too small. We conclude that 
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Fig. 11 Comparison of the linear conductance Fi §- 12 Imaginary part of the pseudo fermion self 

G(T) from FRG (dots) with the result obtained via ener gy r (<^ = °) versus of T / T k (note the loga- 

the NRG |[29l rithmic scale). Inset: T at T/Tpc < lona linear 

scale. 
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Fig. 13 Conductance G(V) at T = from FRG. 
The lines indicate the half width at half maximum 
which amounts to AV = 1.93Tj^. 
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Fig. 14 Imaginary part of the pseudo fermion self 
energy F(uj) = versus V/T^ on a logarithmic 
scale. Inset: T(0) at V/T^ < lona linear scale. 



higher order terms in the (3 -function are not sufficient to cure the problem with T . In addition, the 
frequency dependence of both, g and T is essential in reaching the Fermi liquid regime, as we now sketch. 

First we note that unitarity does not require g(uj) to be bounded at all frequencies. It appears possible to 
have g(u> = 0) growing large or even diverging at strong coupling, while the conductance G remains finite. 
This is because G as given by equ. (f35T > is given by integrating the product of two g 's over energy. The 
energy integration removes any singular behavior of g at u> — provided the singularity is not too strong. 

Secondly, a sufficiently strong relaxation rate T might suppress the unphysical growth of g . As dis- 
cussed in detail in section 3.3 the RG equations in the approximation of taking all energies on shell do not 
generate a sufficiently strong T(uj = 0) at the scale A ~ Tj( , where it would be needed to cut off the 
growth of g , for the reasons discussed above. On the other hand, r at A = turned out to grow too large at 
T <C Tk ■ This deficient behavior may also be remedied by taking the energy dependence of T and g into 
account. Qualitatively one expects that the smearing of the sharp singular structures in g(u>) and T(u>) by 
integrations over u> within the width T specified by the pseudofermion spectral functions will be equivalent 
to replacing g(to — 0) and T(uj = 0) by g(u> ~ T) and r(o; ~ T) , which are both finite quantities of 0(1) 
and 0(Tk), respectively, not only at the end of the RG process (at A = ), but already at A w Tk ■ 
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4 Conclusion and outlook 

With the advent of nanophysics in the mid 1990s it became possible to fabricate nanostmctures with the 
properties of single atoms. In contrast to real atoms these artificial atoms may be contacted by metallic 
leads, opening a new way of probing their properties in transport measurements. This development led to 
a revival of Kondo physics, a subject that had emerged in studies of dilute magnetic alloys. Even though 
the Kondo effect had been well understood by the mid 1970s, the new development required theories 
describing the Kondo effect out of equilibrium. Most of the theoretical methods created to describe the 
Kondo effect in equilibrium do not work or work only in a limited way for systems out of equilibrium. A 
good starting point is perturbation theory in the Keldysh formulation. However, due to the infrared singular 
behavior of the Kondo problem, simple perturbation theory is not sufficient. In equilibrium the powerful 
method of the renormalization group has been proposed and successfully implemented first by Anderson 
in the perturbative regime and later by Wilson, all the way to strong coupling. 

Anderson's perturbative treatment has been generalized to non-equilibrium in iTTOl [TP , after recog- 
nizing two important new features. First, in nonequilibrium the coupling constant has to be generalized 
to an energy dependent coupling function. Second, the processes limiting the quantum coherence of the 
allimportant spin flip processes that are responsible for the Kondo effect in the first place, become more im- 
portant out of equilibrium and may stop the RG flow to strong coupling inside the weak coupling regime. 
The latter lead to a decoherence rate, which in iflOl ITT1 was added ad hoc to the RG equations. In the 
present paper the decoherence rate is treated from the beginning on equal footing with the coupling func- 
tion. It is known |[30l that the decoherence rate is given by the physical spin relaxation rate. In technical 
terms, the spin relaxation rate has contributions from the self energy as well as from vertex corrections. 
Within the standard quantum field theoretical formalism used in this paper, it is straightforward to incor- 
porate the self energy contribution. To include the vertex corrections would require to take into account 
three particle processes as well. The latter are beyond the scope of what has been done in the present 
paper. There exist different formulations based on the Liouville operator formalism 11331 or the flow equa- 
tion method J32), which allow to take into account vertex corrections in a systematic way. However, these 
methods have their own drawbacks. The main achievement of the present paper is therefore the principal 
demonstration of the appearance and the actual calculation of the effect of the decoherence rate on the RG 
flow for the Kondo model. 

It is tempting to explore the importance of the decoherence processes in cutting off the RG flow to 
strong coupling in equilibrium, where it is strongest. From the known behavior of the spin relaxation rate 
as a function of temperature one would expect that at temperatures at and below the Kondo temperature 
the RG flow is already cut off and the coupling constant would tend to a temperature independent value 
with T 2 -corrections at low T . We found something different: apparently the coupling function g(uj = 0) 
grows to values much larger than unity, without, however, violating the unitarity limit. A careful solution 
of the frequency dependent RG equations is a subject for future work. If the coupling grows to large values 
one should expect higher order processes to become important. We have added higher order terms to the 
RG equation, in such a way that g(u> = 0) tends to a finite value of order unity, recovering the unitarity 
limit of the conductance in its simplest form. The resulting conductance curve is in good agreement with 
NRG results known for the equilibrium case, and allows to explore the behavior at finite bias voltage. 
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